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Abstract 

We establish new relations which connect Euclidean sonar transforms 
(integrals taken over spheres with centers in a hyperplane) with classical 
Radon transforms. The relations, stated as operator identities, allow us 
to reduce the inversion of sonar transforms to classical Radon inversion. 

1 Introduction 

As we aim to relate sonar transforms with Radon transforms, we must begin 
by recalling key definitions which in turn requires us first to fix some notation. 
E'l := R"^i X (0, oo) denotes the upper half space of W\ Points in will be 
written as {x, y) with x G ]R"~^ and y > 0. We write \x\ for the Euclidean vector 
norm of x G R"^^, dt for the Euclidean volume element on R""-'^. {^+) 
(resp. C°° (R")) denotes the set of smooth compactly supported functions 
(resp. smooth functions) on R" supported in R" . S""""^ will denote the sphere 
in R" centered at x with radius y (empty if y < 0) carrying area measure dS. 
We now define the sonar transform S : {^+) ^ (^+) follows. 

Definition 1.1. Given / e {W]^), 



S[f]{x,y) := 




The centerset variable x parameterizes the centerset R" ^ x {0}. On occasion 
we call y the radial variable. 



More explicitly, 



S[f]{x,y)=[ fL + t,V^^^) (1) 

J\t\<y ^ ' vr - \A 

The sonar data generally does not have compact support. However the 
restriction of S\f\ to any hyperplane parallel with the centerset is compactly 
supported which justifies various compositions of transforms below. 

Remark 1.2. While we restrict <S to {^+) foi" the sake of our subsequent 
derivations, Definition 1.1 makes sense for locally integrable functions /. 

Courant and Hilbert initiated the study of S in Methods of Mathematical 
Physics, Volume II [2], where they established its injcctivity on the space of 
continuous functions and used the result to investigate hyperbolic partial differ- 
ential equations. While [2] terms the mapping S "integrals over spheres centered 
in the plane," our more efficient "sonar" terminology follows recent applications 
of <S to marine tomography as in work of Louis and Quinto [8] , where the oper- 
ator S (in dimension three) models naval sonar data. 

Operator S has other practical uses. As Cheney [1] explains, in dimension 
two S describes synthetic aperture radar. As S and its generalizations abstract 
the behavior of reflected waves (echoes) whether acoustic, electromagnetic, or 
mechanical, they play a central role in reflective tomography, including marine 
tomography and radar theory. 

For the sake of recalling the classical Radon transform TZ, we will write 
the set of all hyperplanes in a real vector space V as ViV), which carries the 
structure of a smooth manifold. We now define TZy ■ {V) i— > ('P(^)) as 
follows. 

Definition 1.3. 

T^vUKP) := J^f{x)dm,, 

where dm^ denotes planar surface measure on P e ViV). 

For V = M" we write the Radon transform simply as TZ whereas for V = 
]^n-i ^ |Q| (- denote the Radon transform TZ and call it the centerset 

Radon transform. Note that both the sonar and Radon transforms reduce to 
the identity operator when n = I. 

The operators <S and TZ appear quite different conceptually. For example, 
Helgason [6] shows that the Radon transform has a very special group-theoretic 
structure which leads to the following inversion formula in M" (taken from [6], 
P-15): 

Here the superscript * indicates the adjoint operator while A stands for the 

Laplacian. In odd dimensions, (— A)^~ is a differential operator; in even di- 
mensions, the fractional power of the negative Laplacian takes the form of a 
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pseudo differential operator that should be interpreted in terms of Riesz poten- 
tials. 

In sharp contrast, the composition S* oS does not exist, as S[f] may lack 
compact support even when / has one. And, despite the sonar transform's 
sizable symmetry group (on R" , a semidirect product 0(n — 1) ex M""^ of 
the orthogonal group with translations), at present we lack a group-theoretic 
interpretation of S. 

In light of these differences, unexpected close relations between these two 
operators carry intrinsic interest. Denisjuk [3] found the first such relation: 

Theorem 1.4 (Denisjuk, 1999). B^ ;^ denotes the unit ball in M". There 

exists a mapping (j) : B" R" (related to stereographic projection) and a 
certain non-negative weight cr : B" ^ — > R such that for any smooth compactly 
supported function f G (^+) > 

S[f] =n[a- if o <!>)]. (3) 

Using Equation (3), Denisjuk expressed the inverse as a pull-back of the 
inverse Radon transform Ti~^ and established Plancherel identities for sonar. 
Theorem 1.4 was later used in [9] by Palamodov to pull back various microlocal 
estimates and perform A-type reconstruction on integrals over arcs. As both [3] 
and [9] demonstrate, sonar-Radon relations of type (3) can be effectively used 
to translate any Radon result into a corresponding sonar statement. 

Denisjuk connects the sonar transform of a given function to the Radon 
transform of a different function with, a priori, a different support. Since 
Rl C R" and thus (R^) C (R") , it makes sense to speak of the Radon 
transform and sonar transform of one and the same function / (as long as / has 
support in R" ). This article addresses the following natural question: 

How can one pass directly from the sonar transform of a function to 
its Radon transform? 

A simplified version of our main result says that for almost all planes P, we 
can compute TZ[f]{P) as 

{WoA^fyov^ onoS)[f]iP). 

for certain explicitly described operators W, Ai/y, and T>ri-2 each of which has 
a natural geometric meaning. 

2 Main Result 

We aim to compute the function TZ[f] explicitly from the function S[f] , the sonar 
data associated to /. In order to make our calculations as explicit as possible, 
we need a suitable parameterization of P(R"). Our formula for calculating 
TZ[f]{P), as it turns out, breaks into various cases depending on the geometry of 
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the hyperplane P relative to the centerset of the sonar transform. Accordingly, 
we write 



a disjoint union, with 

h the set of planes parallel to M""^; 

V the set of planes perpendieiilar to M"^^: and 

s the set of all other planes (the slanted planes). 

We shall write Ti^'^^[f], 7i^'"^[f], and 7^*^*')[/] for the corresponding restrictions 
of TZ[f] to h, V and s. 

To make matters more precise, we begin with a general parameterization 
of P. To a pair {u>,p), w € P — 0; "^^J associate the hyperplane 

P = {x e W^\u) ■ X = p}; here • denotes the standard inner product in M". 
P almost determines {uj^p)\ only if p = 0, w may vary by a sign. 

Now we adapt this framework to take account of the centerset. We write 
uj = e x M, and similarly x = {x',Xn)- The equation oj ■ x = p 

now takes the form w' ■ x' + cOn Xn = P- We can now distinguish three cases, as 
above: 

{h) co' vanishes; 

(v) uin vanishes; 

(s) all others. 

Start with (h). Since iv e we must h ave LOn = ±1. Dividing through by 

oJn, the defining equation of the plane has the form Xn = y (for some appropriate 
y). The variable y can now parameterize {h). If / € (^+) then 



Thus we integrate over the centerset variable. 

Now turn to {v). A plane in v has defining equation iv' ■ x' = p. The inter- 
section of a vertical plane with the centerset determines it, so parameterizing v 
amounts to parameterizing the set of all hyperplanes in the centerset. As before, 
P almost determines (w',p); only if p = 0, iv' may vary by a sign. Accordingly, 
we shall write 



ViW) = h[Jv[Js, 




(4) 




(5) 



where w = (w', 0). Observe that 




(6) 
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(which we view as an equality between functions of y), so we may also write 
Equation (5) as 

/•oo 

n^''\f]{u:,p)= / n[f{;y)]{uj,p)dy. 
Jo 

The final case (s) comprises a dense open subset of T'(R"), and thus makes 
the most substantial contribution to our union (provided n > 1). 

A plane in (s) has defining equation u)' -x' +WnXn = P with neither lo' nor a;„ 
vanishing. We can scale this equation so as to normalize oj' and simultaneously 
render the coefiicient of Xn negative. Thus we can unambiguously choose a 
defining equation for the same plane with 

(uj',p, -w„) e §^7^ X M X 1R+ 

Now the pair {cu',p) by itself determines the intersection of P with the centerset, 
so ojn controls the angle between P and the centerset. More explicitly, write 
u}n = — cot p with (3 e (0, 7r/2). Then the defining equation of the plane has the 
form 

Oj' ■ x' — cot (iXn = P- 

Intersecting with the parallel translate of the centerset where a;„ = s sin f3 gives 

uj' ■ x' — p + s COS l3. 

Henceforth we parameterize (s) by (w',p, (3 now directly represents the angle 
between P and the centerset. 

Integration of / G (^+) ^^^^ P can be split into integration over P fl 
followed by integration over (3. Explicitly, 

n('\f]{J,p,p)= n[ f {; s sin (3)] {uj,p + s cos f3)ds. (7) 
Jo 

Remark 2.1. We avoid including (v) in (s) as a special case (3 = 7r/2 both to get 
a good parameterization of (s) and because the cases require separate treatment 
below. 

Theorem 2.2 (Sonar-Radon relations in R"). The sonar transform S[f] 
determines TZ[f] by means of the following operator identities: 

= Dn-i o TZ'^''^ o S, (8) 

n^"^ = CoS, (9) 
= WoAi/yoVn^olloS. (10) 

Here TZ stands for the centerset Radon transform; Aa{y) stands for the weighted 

Radon transform, (38) from Definition 7.1 in Section 7; I?^ and W denote frac- 
tional operators defined in Section 3 by (18) and (20), respectively; L represents 
an infinite limit defined in Section 7 by Equation (33). 
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Wc organize the proof as follows. Sections 3 and 4 contain necessary an- 
alytical tools: Section 3 details the fractional operators and W (and their 
inverses); Section 4 collates identities for spherical integrals of plane waves from 
F. John's classic Plane Waves and Spherieal Means [7], for use in Section 8. 
Section 5 and 6 treat TZ^'^^ and TZ^'"\ respectively. Sections 7 and 8 treat TZ^^^ 
and thus complete the proof of our sonar-Radon relations: Section 7 motivates 
the choice of the weight a = 1/y for operator Ai/y and establishes results in 
dimension two; Section 8 generalizes these results to higher dimensions. Section 
9 analyzes the main result and offers closing remarks. 



3 Fractional Calculus 

We recall notions from fractional calculus, especially regarding operators 7)^ 
and yy appearing in our sonar-Radon relations. The standing assumption that 
all fractional operators act on the last variable of smooth compactly supported 
functions will avoid those various delicate issues discussed at length by Samko et 
al in [10]. Compact support obviates potential divergence; smoothness ensures 
commutativity of operators. 

As we choose to view the fractional integrals 1^ as the fundamental fractional 
operators, we develop all other fractional operators out of these. 

Definition 3.1. For i/ > we define X^, a type of fractional integral, by 

U9]{y) ='^f j\y'- ais)ds. (11) 

For convenience, we also set To = id ., the identity operator. 

According to Lemma 3.2 below, the set {Ii^y^^Q of fractional integral oper- 
ators forms a monoid (semigroup with identity) under composition. 

Lemma 3.2. The fractional integrals in Definition 3.1 satisfy the composition 
law 

1tJ.°1v =Iix+v, (12) 

which holds for all /x, > 0. 

Proof. Equation (12) certainly holds if either = or = on account of the 
convention Iq =id. So assume /x, > 0. 

Using Equation (11), we express the composition {I^oIi,)[g\[y) as an iterated 
integral 



- [y -s) s [s -t) 
Jo Uo 



g{t) dt 



ds. (13) 



Changing the order of integration in (13) yields the following expression for 

{I^oI,)[g\{y): 

^^^y[m]mdt. (14) 
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where 



Substituting = {y^ - t^) P + t^, gives 



2s ds. 



(15) 



with the constant 



C = [ {1 - p)"-^ p"-^ 
Jo 



dp 



taking the form of Euler's integral of the first kind (see [4] p. 948 # 8.380.1) with 
value given by 



C = 



Thus 



2 7^^*+^ ry 



r(M)r(;^) / 2_ .2\/^+'^-i 



g{t)dt. (16) 



Cancelhng four gamma terms, we obtain 



as desired. 



□ 



The justification for the terminology fractional integrals for the lu rests on 
the semigroup property and the observation that 



rv 

^ib](y) = 2 7ry / g{s)ds, 
Jo 



a scaled antiderivative. Now Ii admits a left inverse in the form 



'Di[g]{y) :- 



1 d 



9{y) 



y 



2iT dy 

So, symbolically, we have I?i o Xi = Xq = id, and, more generally. 
Lemma 3.3. For all ly > 1, Vi oT,y = 

Proof. Using Lemma 3.2 and the associativity of operator composition 
©1 o = Di o (Ji o J^_i) = (Di o Ii) o = J^_l. 



(17) 



□ 



We will now use Lemmas 3.2 and 3.3 to construct fractional derivatives of 
arbitrary order. 



Definition 3.4. For u > the fractional derivative V^, is defined in terms of 
(11) and (17) as a mapping 

V^ = V[''^ oI^,^_,, (18) 

where is the smaUcst integer greater than or equal ly. For consistency with 
Definition 3.1, wc define Vq = id . 

As an immediate consequence of Lemmas 3.2 and 3.3, we have the following 
corollary. 

Corollary 3.5. For u > the fractional derivative is the inverse of the 
fractional integral!,,, i.e.: "D^ oX^ = id. 

In Section 7, we will encounter a fractional operator V and require its inverse 
W to deduce Equation (10) in Theorem 2.2. 

Definition 3.6. Set 

^^am--= f ^. , mde. (19) 

^0 v sin p — sm 6 
From the definitions of Ii and V 

2 

Ji[5](sin/3)=V[cos/35(sin/3)](/3). 

We now cast this identity of functions as an identity of operators. Define 
Q[(/](/3) := g{sm.f3) and K[g\{f3) := cos/3 • g{sinj3). Then the identity above 
says 

Qo Ji = Vo/C. 

2 

Thus 

W := V"^ =/Cor)i oQ-i. 

2 

Lemma 3.7. 

m9m--=-j-J ^=^^==9{e)de (20) 

dp Jo V sin^ j3 - sin^ 6 

Proof. We apply the composition of operators /C o Vi o to a function g. 

First, 

Q~^[9]{s) = 5f(arcsins). 
Next we use Equation (18) to write 

(Pi o Q-')[gm = (Pr o J. o Q-')[gm 



J_ A. 

2^ dp 

1 d g((arcsins) , 
' -ds 



-.2/3/ ii^^ds 



P Jo 



2 



ld_r 
n dp Jo 



.2 
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Finally, we apply the operator K, which replaces (3 with sin/3 and multiplies 
the result by cos/3, to get 

r.. ^ ^ ur i/^x ^1 d r'"'^ 5(arcsins) , 

(KoVi o Q-^)[g\m = cos 13--— — ., ' ds 

whereupon the substitution s = sin ^ yields the statement of the lemma. □ 



4 Plane waves and spherical means 

The identities for integrals over spheres and balls in M" collected here, combined 
with the formulae from Section 3, form the crux of the derivations presented in 
Sections 5 and 8. In particular, we state the Co-area Formula, following [5], and 
develop some of its consequences. 

Theorem 4.1 (Co-area formula). Let u : i-^- M 6e Lipschitz continuous 

and assume that for almost every r S R the level set 

{x e M" I u{x) = r} 

is a smooth, (n — 1) -dimensional surface in M". Suppose that also g : M" i— > R 
is continuous and locally integrable. Then 



J : 



poo 




grdidu\dx = 


f gdS 


J — oo 


J {u=r} 



dr. 



where dS denotes surface measure on the level set {u = r}. 



(21) 



By setting u = \x\ in Theorem 4.1, one obtains a standard identity for 
converting integrals over balls into integrals over spheres in R" which we state 
in Lemma 4.2. 

Lemma 4.2 (Polar Coordinates). Let g : B"^ i— > R be a continuous function 
on a ball of radius r in Euclidean space. Then 



/ g{x) dx= / f{x) dSx 

J\x\<r Jo JxeS"p^ 



dp, 



(22) 



where dSx denotes surface measure on the sphere S"~^ of radius p. 

Differentiating Equation (22) gives rise to the following. 
Corollary 4.3. Let 5 : ^ i-^- R be a continuous function. Then 



d 
dr 



g{x) dx 



g{x) dSx 



(23) 



holds for all < r < R. 
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Often it is convenient to replace integration over a sphere of radius r with 
integration over a unit sphere. As we remark below, this can be accomplished 
through a simple substitution. 

Remark 4.4. Let 5 : S"^ i— > R be a continuous function on a sphere of radius r. 
Then 

/ g{x)dS^ = r^-^ I 9{re)dne, (24) 

where dQ.e is the surface measure on a unit sphere. 

Using Remark 4.4, we recast Equation (22) in the form we shall find most 
useful: 



/ g{x)dx= [ s"-^ / g{s0)dne 
J\x\<r Jo -^fes",!^ 



ds. (25) 



Throughout the rest of this section g denotes a continuous function of a 
scalar variable. Fix v e M". Following F. John in [7], we call G{x) := g{v ■ x) a 
plane wave with normal v; such a function is constant on planes perpendicular 
to V. 

We follow F. John in [7] to reduce integrals of plane waves over spheres and 
balls to single-dimensional integrals. On the plane v ■ x = \v\p the plane wave 
g{v ■ x) has constant value (/(|f|p); the intersection of that plane with the ball 
of radius r forms a ball of radius ^/r^ — p^. Thus 



|x|<r 



g(v-x)dx= \ No\{W-' )g{\v\p)dp 

26 

|gn-2| .+r ^ ^ 

= ^'^^ J {r^-p^) ' 9{\v\p)dp. 

whore denotes the total surface measure of a unit sphere in M"^"'^. Dif- 

ferentiation of Equation (26) with respect to r followed by evaluation at r = 1, 
leads to the following fundamental identity (c.f. [7], p. 8). 

Theorem 4.5. Let w e M" he a fixed vector and let g : R M. be a continuous 
function. Then 



Jees" 



g{ve)dne = \S:;,^\ r\l-p')"~^' g{\v\p)dp. (27) 



We state two consequences of Theorem 4.5. If v has unit length, then Equa- 
tion (27) becomes an identity for "spherical plane waves" : 

Corollary 4.6. Let lj g §"7^^ fixed point on a unit sphere in K" and let 
3 : R i-^- R 6e a continuous function. Then 

g{u-e)dile = Kf\ r'il-p^)"^ g{p)dp. (28) 
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Alternatively, on setting g = I, Theorem 4.5 gives a recursion for the total 
measure of a unit sphere (c.f. [7], p. 9). From this recursion follows the well- 
known surface area formula: 

\Ki'\ = ^- (29) 

We shall use this formula to connect our fractional integrals with geometric 
transforms. 



5 Integrals over horizontal planes 



We now prove Equation (4) of our Main Theorem; 

Te^'') = Vr^ o 'RS^'^ o S. 

2 

Since inverts the fractional integral 1^, it suffices to show that 

n'-^Ks = i„-i on^''\ (30) 

2 

We compare (7^('') o5)[/] with (I^-i oTZ^^'>)[f] for / G (M!^) . Directly from 
the definitions of the operators S and TZ^'^^ (Equations (1),(4)) 

idt 

/xeM"-i lJ\t\<y ^ ' 
which, upon interchanging the order of integration, equals 



(7^(^)o5)[/](y)= / / fU + t,Vy^^) 

JxeM"-i J\t\<v ^ ' \ - 



dx 



J\t\<y LJxGR"-! ^ ^ 



dx 



ydt 



We recognize the inner integral as Ti-^^^ [f]{\^ — and deduce 

ydt 



(31) 



The right-hand side of Equation (31) is an integral over a ball in R" ^ whose 
integrand 

««[/l(v/7^)^^ 

is a radial function of t. Therefore, using Equation (25) from Section 4, we 
rewrite Equation (31) in the form 



Jo 



■ dr. 



which, after substituting r = \/y^ — s'^ becomes 

\Kf\y f {v' - n(^\f]{s) ds = (i^ o 7^("))[/] 

by Definition 3.1 in Section 3 and (29). 



(32) 
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6 Integrals over vertical planes 

Figure 1 suggests viewing a vertical hyperplane as a limiting case of expanding 
tangent spheres with a fixed point of tangency located on the centerset. 
Accordingly, 

7^('')[5](c.,p)= lim S[g]iu;s,\s-p\). 

I s|— >oc 

If we make the definition 

C[f]{Lu,p):^ ]im f{L,s,\p-s\). (33) 

\s\ — 'OC 

then TZ^^^ = Co S, as desired. 



Figure 1: Vertical rays as limits of arcs 



7 Integrals over slanted lines 

In M^, our desired sonar- Radon relation Equation (10) reduces to: 

n^''^ ^WoAi/yoUoS. (34) 

(As T>ri-2 = T>o = id, no fractional derivative appears.) Below, this formula 
emerges as the foundation for the general case. 

In dimension two, the centerset has dimension one. As hyperplanes in di- 
mension one coincide with points, just in this section we will encode them as 
such (rather than as pairs of a unit vector and a magnitude) . With this encod- 
ing the centerset Radon transform TZ reduces to the identity map, so we must 
prove that 

7^(^) = W o Ai/y o S. (35) 
Applying V, the inverse of W to both sides yields the equivalent statement 

V o 7^(^) = Ai/y o S (36) 

for which we will aim. 

Dimension two affords us a simple formula for the sonar transform: 

S[f]{x,y)=l f{x + y cos (t),y sin (j))yd(l}. (37) 
Jo 
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We now furnish the definition of Aa — a type of weighted Radon transform on 
Rl. 

Definition 7.1. Fix any set T. Consider a function g on T x such that 
g{uj,-) e {^+) fo'' 6ach uj in T. For /3 e (0, 7r/2) and non-negative weight 
tr : i-^ R_|_, define a weighted Radon transform by 

/•OC 

y^crbKi^jP, /3) = / .9 + s, s sin/3) cr(s) ds. (38) 
Jo 

This section has T a singleton and we thus suppress the variable uj. 

Consider the composition Aa-°S for a general weight a. S[f] (p, y) means the 
integral of / over over a radius y circle centered at p on the a;-axis. By definition, 
the operator Aa integrates functions along rays with slope sin (3. This makes 
./lcr[5[/] (3) a weighted integral of integrals of / over a family of circles, as in 
Figure 2. 



Figure 2: Semicircles tangential to a ray 



One also sees from the figure that arcs with apexes on slanted rays sweep 
infinite wedges. If the apexes lie on a ray with slope sin [3 then the corresponding 
wedge has slope tan/3. Therefore ^cr[5[/] (3) can be expressed as an integral 
over an infinite wedge with vertex at p on the x-axis and angular measure /?. 
Explicitly, by means of Equations (37) and (38), 



Aa[S[J]]{p,l3) 

f {p + t{l + su\(3 cos (/)) , t sin (3 sin c 



t sin (3 a{t) dt 



(39) 



We now make a change of variables designed to simplify the argument of / in 
Equation (39). Define * : R^ -> R^ by 





* 


p + i (1 + sin/3 cos 


t 


1 — > 


t sin/3 sint/f 
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Observe that 5" sends a line segment connecting (0, t) and (tt, t) to the semicircle 
centered at {p + t, 0) with radius t sin/3. 
On each of the two semi-infinite strips 

(O, I + /?) X (0, oo) and (| + /3, tt) x (0, oo) 

\E' acts as a diffeomorphism to the infinite wedge 

{(x,y) € I p < X, <y < tan/? (a; - p)} 

shown on Figure 2. In terms of \E', the double integral (39) over a wedge can be 
written as a sum of two integrals over infinite strips 



/I/ 

Jo Jo 



(/ovI/)(^,f)d<^ 



i+0 



t sin (3 a {t) dt 



t sin (}(j(t) dt. 



(40) 



Introducing polar coordinates (p, 9), p G (0, oo), 6 G (0, (3) in the wedge 

{p + t{l + sin(3 cos(f)),t sin/3 sin^) = (p + p cos6,p sinO) 
gives us the relations: 

t(l -l-sin/3 COS0) = p cos^ (41) 
t sin /3 sin = p sin 6. (42) 

Using (41) and (42), we shall now change (40) into a much more amenable 
expression. 

In order to transform (40), we need to express the old variables {(p,t) in 
terms of the new variables {6, p) and find the corresponding Jacobian factors: 
one for each integral in (40). From the algebraic point of view, it is easier to find 
(j). Divide Equation (42) by Equation (41): this eliminates variables p and t. 
Next use trigonometric identities to solve the resulting relation between angles 
as follows: 

Eq. (42)^ sin 13 sin 4> _ ^^^^^ 
Eq. (41) 1 -l-sin/3 coscf) ~ 

sin /3 sin (j) cos ^ = sin ^ + sin /3 cos (j) sin 6 

sinO 



sin 6 cos 9 — cos 6 sin 6 = 



sin/3 



sin(^ -9) = 



sme 



sin/3 

we get two solutions: 

' sin9 



(43) 



' + sin 



sin/3 



TT — sm 



sin/3 



> . 
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With these expressions for the angular variable 0, we may now find the corre- 
sponding values of t as outlined in Equation (44) below: 

Eq. (41) X cos 6* + Eq. (42) x sin 6 

=> t cos 9 + t sin /? (cos cf) cos 9 + sni(f) sin 9) = p 
t cos 6 + t sinp cos{(f) — 9) = p 

^t= 

cos 9 + sin (3 cos((/) — 9) 

using the identity following from Equation (43): 

sin /3 cos(0 - 61) = ± sin /? ^ 1 - sm^ {(j) - 9) =\ (44) 

^± sin/3 -^1 - (sin^ 6/ sin^ (3) = i-^/sin^ /3 - sin^ 6 j 

we get two solutions: 

P 

cos 9 + \/sin^ 3 - sin^ 9 ' 

P 

cos — \/ sin^ — sin^ 6 

Remark 7.2. The relation between {(p,t) and (^, p) can also be derived geomet- 
rically. 





O A 

Figure 3: Polar coordinates 



Consider semicircles inscribed in a fixed wedge of angular measure /3 as in 
Figure 3. Our old variables (0, t) specify a semicircle and then a point B on it: 
from t = OA we learn the center of the semicircle, and then from the tangency 
also its radius; cj) locates the point B since Z.OAB = n — (p. From the right 
triangle AO AC, we find the radius of the semicircle r = t sin/3. A ray issuing 
from O at angle 9 < (3 with OA will meet the semicircle twice and we take B 
as the second intersection. The Law of Sines applied to triangle AOAB gives: 

sm{A AOB) _sm{Z.OB A) sin 6' _sin(^-6i) 
AB ~ OA ^ t sin/3 ~ t ' 
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which is equivalent to (43). Then the Law of Cosines, in the form 
OA^ + AB'^ -20AAB cos{AOAB) = OB^, 

tells us 

+ sin^ /3 - 2 sin /3 cos(7r -</.)= p^. 
After simplification (44) results. 

As follows from Equation (43) the angular variable (j) docs not depend on p. 
Therefore the 2-by-2 Jacobian matrix is triangular and its determinant is given 

by 

d(h dt 



det 



£2 

de dp 



The values of the partial derivatives -j 
through straightforward differentiation: 



" de ^ dp' 

and ^ for i = 1,2 can be found 



d(l>i 
86 



i-iy 



COS! 



y/sin^ 13 - sin^ 9 + {-ly cost 



^l-(sin2e/sin2/3) ^^^^^ ^sin^ /? - sin^ 

5ti 1 
dp 



cos I 



+ (-l)Vsin^/?- sin^fi'' 



i = l,2. 



whence follows that for both sets of variables ti) i = 1,2 the absolute value 
of the determinant of the Jacobian is given by the same simple expression 

1 



Vsin^/^-sin^fil" 

We conclude that (40) can be written as a single integral of the form 



sin/3 {tia{ti)+t2<j(t2)} 
s/sm^ (3 - sin^ 9 



f{P + P cos 9, p sin 9) dp 



d9, 



(45) 



where the values of t^, i = 1, 2 in the numerator are given by (44). 

Integral (45), representing the composition {Ac,{y) o S) applied to /, becomes 
particularly simple if one sets the weight a = 1/y : 



2 sin/3 



-L 



\J sin^ /3 — sin' 
2 sin/3 



,f{p + p cos 9, p sin 9) dp 



d9 



1/ sin^ /3 — sin^ 9 



/ f{p + pcos9,psm9)dp 
Jo 



(46) 



Recognizing the bracketed integral as TZ^'^^ (via Equation (7) from Section 2) 
and noticing that the outside integral is the fractional operator V (Definition 
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3.6 from Section 3), we get 



{AyyOS)[f]{p,(3)= [ 

Jo 




as desired. 

8 Integrals over slanted planes 

We now prove in all dimensions the sonar-Radon relation (10) first stated in 
Theorem 2.2 and reproduced below: 



After some work, we reduce to the two-dimensional case treated in Section 7. 
EfTcctivcly, our coim^rsion of sonar data into integrals over hyperplanes proceeds 
through an intermediate stage — integrals over cylinders. 

In the usual way, let {lj,p) encode a hyperplane in the centerset of . By 
a cylinder, with radius r with axis we mean any set: 



We encode a cylinder of radius y as a triple {u),p,y) and write C[f]{uj,p,y) for 
the integral of / over the given cylinder. (One naturally views transform C as a 
hybrid of sonar and Radon.) Given S[f], we can find C[f]{LO,p, y) as follows. 



Proof. We shall actually prove the equivalent claim l^iqi oC = TZoS. Combining 
Equation (1) from Section 1 with Equation (6) from Section 2, we obtain an 
iterated integral for {TZ o S)[f]{w,p, y) in the form: 



Interchanging the order of integration, which is possible because / is smooth 
and compactly supported, we get 



7e(«) = W o Ax/y o Vn^ oTloS. 



{{x,y)eRl I {wx-pf+y^=r^}. 



Theorem 8.1. For f G {M.X) , 



C[f]{uj,p,y) = {Vn^o-RoS)[f]{w,p,y) 



(47) 
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where the inside integral is a Radon transform of a shifted function: 
/ f Ix + t, \Jy'^ — dm^ (substitute x-\-t = u) 

J UJ-X=p 

= I fU Vy^ - \t\A dmu 

We conclude that (7?. o 5)[/](a;,p, y) is the following integral over a ball 

J\t\<y ^ ^ 



which, after switching to polar coordinates (Equation (25) from Section 4), 
becomes 



»ra-2 



dr. 



Inside the brackets, we have an integral of a plane wave over a unit sphere. 
Therefore, in Hght of Corollary 4.6, the composition {TZ o S)[f]{(j,p,y) can be 
expressed as the following double integral: 

\K:i'\y / ^ (1-^') = dsdr. (48) 

Jo J-i \/y —r-' 



The mapping 



(s, r) (p + r s, ^/y^ — j 



is a diffeomorphism from the rectangle [—1,1] x [0, y\ into an upper half-disk of 
radius y centered at p. This suggests the following change of variables 



V 



-\-r s=p-\-u, — = 



where v G [0, y] and u < \v\. 

Wc will now transform the integral in (48). Solving for (s,r) in terms of 
{u, v), we find that 



s = 



sj y^ - -V iP' 
and therefore the integrand 



„ri-2 



(1 - S-) 



2\ — 
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in Equation (48) transforms to 

(y2 _ ^2^ ^ (^^^ p^^^ _ ^2^ 

As can be found through routine differentiation, the Jacobian of the mapping 



y2 _ 



{u,v) 



is the following matrix 



(y2_^2_,_„2)3/2 (y2_^2+„2)3/2 



2_,,2+„2-il/2 (j,2_^2_)_„2)l/2 



whose determinant has the absolute value 

V 



Now (48), representing (7^ o 5)[/](a;,p, y), takes the form 



du 



dv. 



From Definition 3.1 of 1,^ (and the surface area formula for spheres found at the 
end of Section 4) we may conclude that 



{n o S) [/] {w, p, y)=ln^lj^\[f](uj,p + u, - u^) 



y du 



But now we recognize the expression in curly brackets asC[/](a;,p, y), as desired. 

□ 

To finish, note that fixing u) determines a parallel family of cylinders. In the 

definition of Aa we now set T equal to the set all possible w, i.e. T = §"7^- 
According to Section 7, the composition W o Ai/yoC yields the two-dimensional 
Radon transform of TZ[f] which is the n-dimensional Radon transform T?,*-*-* [/]. 

9 Conclusion 

The technique which proves the main theorem admits immediate variations, 
if perhaps of only theoretical interest. For the record, we mention two. In the 
sonar transform one could replace the spheres that function as loci of integration 
by other families of loci with similar scaling properties. Alternatively, in the 
Radon transform, one could replace slanted planes by cones whose axes lie in 
the centerset. 
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The authors view the methods in this paper as an expression of a more gen- 
eral philosophy, under development, aimed at providing sonar-Radon relations 
for more general centersets and in more general spaces. The planar centerset 
case deserves an independent treatment now because its rich structure allows 
for results of a particularly explicit form and because of potential for practical 
applications. 
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